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ABSTRACT 

Estimating the black hole mass at the center of galaxies is a fundamental step not only for un- 
derstanding the physics of accretion, but also for the cosmological evolution of galaxies. Recently a 
new method, based solely on X-ray data, was successfully applied to determine the black hole mass in 
Galactic systems. Since X-rays are thought to be produced via Comptonization process both in stellar 
and supermassive black holes, in principle, the same method may be applied to estimate the mass in 
supermassive black holes. In this work we test this hypothesis by performing a systematic analysis 
of a sample of AGNs, whose black hole mass has been already determined via reverberation mapping 
and which possess high quality XMM-Newton archival data. The good agreement obtained between 
the black hole masses derived with this novel scaling technique and the reverberation mapping values 
suggests that this method is robust and works equally well on stellar and supermassive black holes, 
making it a truly scale-independent technique for black hole determination. 
Subject headings: Galaxies: active - Galaxies: nuclei - X-rays: galaxies 



1. INTRODUCTION 

There is now convincing evidence that the most pow- 
erful persistent sources in our Galaxy (Galactic black 
hole systems; hereafter GBHs) and in the universe at 
large (active galactic nuclei; AGNs) are powered by grav- 
itational accretion onto stellar and supermassive black 
holes, respectively. One of the most evident manifesta- 
tions of black hole systems (BHs) is the presence of rela- 
tivistic jets, which may have a huge impact on their en- 
vironment over distances that are well beyond the radius 
of influence of the black hole itself, as demonstrated by 
recent work on the e ffects of powerful .jets at the center of 
galax y clusters (e.g. iMcNamara et al.ll2000t iFabian et all 
l2003l ). In addition, and perhaps more importantly, the 
tight correlations found between the BH mass (Mbh) 
and several galaxy parameters such as the velocity dis- 
persion or the mass of the bu lge (e.g. iGebhardt et al.l 
l2000tlFbrrarese &: Merritti2000( ) clearly demonstrate that 
black hole growth and the build-up process of galaxies 
are closely related and therefore black holes are essential 
ingredients in the evolution of galaxies. 

The estimate of the mass of supermassive black holes 
at the center of galaxies is the most crucial parame- 
ter needed to understand the formation and evolution 
processes in galaxies as well as the central engine in 
AGNs. Different methods have been developed to mea- 



sure A/bh in AGNs, depending on their characteristics 
(e.g., degree of activity, variability, distance). For ex- 
ample, in the case of weakly active or quiescent galax- 
ies, Mbh can be determined by directly modeling the 
dynamics of gas or stars in the vicinity of the black 
hole (e.g. iKormendv fc Richstonel[l995l : iMagorrian et al.l 
|1998() . On the other hand, for highly active broad-lined 
galaxies that show strong variability, the most widely 
used way to estimate Mbh relie s upon a technic|ue kn own 
as reverberation map ping (e.g., [P eterson et al. 2004' and 
references therein; see lCzernv fc N ikolaiuk 2010 for a re- 
cent review of the most commonly used methods). 

Mass estimates obtained by these methods have pro- 
vided very important results for our understanding of BH 
formation and galaxy evolution. However, all methods 
have specific limitations. For example, methods based 
on the modeling of gas and stellar dynamics require the 
sphere of influence of the BH to be resolved. There- 
fore, these techniques can successfully be applied only to 
nearby galaxies with relatively massive BHs and where 
the galaxy's optical emission is not substantially affected 
by the black hole activity. On the other hand, the rever- 
beration method is time intensive, and cannot be applied 
to very luminous sources, whose variability is typically 
characterized by small- amplitude flux changes occurring 
on very long timescales. Furthermore, this technique re- 
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quires the presence of a detected broad line region (BLR), 
whose nature and geometry is still poorly known. Similar 
limitations affect most of the secondary methods, which 
rely on some empirical relationship between Mbh and 
differ ent properties of the host galaxy (see iVestergaard. 
l2009l for a recent review). Additionally, bright type 1 
AGNs (i.e. the AGNs characterized by the presence of a 
BLR) comprise only a small minority of the total galaxy 
population, severely limiting the number of galaxies that 
can be probed with this technique. Indeed, even amongst 
the population of known optically-identified AGNs, type 
2 AGNs (i.e., the AGNs without visible BLR) appear 
to be about 4 tinies mo re numerous than type 1 AGN 
(|Maiolino fc RiekdHool . 

Exploring the SMBH mass function in as diverse as 
possible galactic environments is essential to understand 
the evolutionary history of BHs and their connection to 
their host galaxies. To address this serious deficiency, 
it is critical to explore alternative methods to determine 
the BH mass that are not dependent on measurements 
at optical wavelengths. 

X-ray observations are one of the most effective means 
to investigate the properties of AGNs for several reasons. 
First, the X-ray emission is one of the defining properties 
of AGNs: while optical lines generally attributed to AGN 
activity may be obscured, suppressed, or undetected for 
several different reasons (e.g., low signal-to-noise ratio 
S/N, excessive redshift, galaxy or star light contamina- 
tion), the X-ray emission appears to be ubiquitous in 
AGNs. Second, the X-rays are produced and reprocessed 
in the inner, hottest nuclear regions of the source. There- 
fore, unlike the optical lines from the broad and narrow 
line regions that are produced by the reprocessing of the 
primary emission, the X-ray emission directly traces the 
black hole activity. Finally, the penetrating power of 
(hard) X-rays allows them to carry information from the 
inner core regions without being substantially affected 
by absorption. 

At first order, from the spectral point of view, the 
AGN X-ray emission in the 2-10 keV energy range is ad- 
equately described by a simple power-law model. Impor- 
tantly, the photon index F appears to correlate directly 
with the accretion rate, suggesting that X-ray spectral 
properties provide indirect information on the accretion 
state in AGNs. While there is growing evidence of a 
positive correlatio n between F and E ddington ra tio for 
brigh t AGNs (e.g. lShemmer et al1l2006; Papadakis et al.l 
|2009| ) , there is also suggestive evidence of a possible anti- 
correlation fo£^I^_Sifsteins accretin g at low Eddington 
ratio (e.g. IConstantin et al.l 2009; W u fc Gul[2008l ). Al- 
though several details are still unknown, it is now widely 
accepted that the bulk of the X-ray emission both in 
AGNs and GBHs is produced by the Comptonization 
process, i.e. the repeated inverse Compton scattering, 
where seed photons likely produced via thermal emission 
by an accretion disk are up-scattered by highly energetic 
electrons. Therefore, X-rays represent one of the most 
effective tools to probe the properties of black holes sys- 
tems and may provide independent constraints on their 
mass. 

In this work we investigate whether a novel method 
to estimate Mbh, which is based solely on X-ray spec- 
tral data and was successfully applied to GBHs by Sha- 
poshnikov & Titarchuk (2009; hereafter ST09), can be 



extended to AGNs. To this end, we utilize a sample 
of AGNs with A/bh already well constrained via rever- 
beration mapping and with good-quality X-ray data. In 
Section 2 we summarize the main characteristics of this 
method. The AGN sample and the data reduction are 
described in Section 3. The results are presented in Sec- 
tion 4 and discussed in Section 5, where we draw our 
main conclusions. 

2. A NEW METHOD TO ESTIMATE Mbh 

Recently IShaposhnikov &: Titarchufl ()2009() carried 
out an extensive and systematic analysis of the tempo- 
ral and spectral X-ray data from 17 spectral transition 
episodes in 8 different GBHs. Their main findings can 
be summarized as follows: (1) Two positive correlations 
are found: the first one involving temporal and spectral 
properties and specifically relating the quasi periodic os- 
cillation (QPO) frequency and the photon index F; the 
second one relating two spectral parameters of a specific 
Comptonization model described below: the normaliza- 
tion A'bmc and F. (2) Both spectral evolution trends 
(QPO - F and A^bmc - E) can be adequately parametrized 
by analytical functions, which are similar for the different 
GBHs. (3) The self-similarity of these pair of trends in 
different GBHs allows the estimate of Mbh and distance 
for these binary systems. 

Since the Comptonization process producing the X-ray 
emission is widely thought to work in the same way in 
stellar and supermassive black holes, in principle, the 
same method may be applied to estimate the mass of su- 
permassive black holes. Clearly, it is not possible to use 
the QPO - F diagram, since no firm detections of QPOs 
have been reported in AGNs (with the exception of the 
NLSl RE J1034+396 in one occasion; iGierlihski et al.l 
l2008f) . But this does not represent a problem for AGNs, 
since their distance is independently and robustly deter- 
mined from redshift measurements or from variable star 
methods for objects of the Local Group. Hereafter, we 
will therefore focus solely on the iVBMC - E diagram. 

In the following, we first describe the general details of 
the Comptonization model used to fit the X-ray spectra 
and the basic characteristics of the scaling method that 
allows one to estimate the mass for any BH system by 
simply scaling the Mbh value from a suitable Galactic 
reference source. 

2.1. The BMC model 

Although the geometry and nature of the Comptoniz- 
ing region are still matter of debate because of spectral 
degeneracy, for our purpose (i.e. estimating black hole 
masses), the specific physical condition of the Comp- 
tonizing medium (i.e. energy and spatial distribution 
of the upscattering electrons and their optical depth) are 
negligible, provided that the Comptonization process is 
characterized in a simple, general, yet physically correct 
way. These criteria are met by the so called Bulk Motion 
Comptonization model (BMC), which is a generic Comp- 
tonization model able to describe equally well the ther- 
mal Comptonization (i.e. the inverse Compton scattering 
produced by electrons with a Maxwellian energy distri- 
bution) and the bulk motion Comptonization (where the 
seed photons are scattered off electrons with bulk rela- 
tivistic motion), although it was historically developed 
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to describe the Comptonization of thermal seed photon s 
by a relativistic converging flow ()Titarchuk et al.iri997[ ). 

The BMC model is characterized by 4 free parame- 
ters: 1) the temperature of the thermal seed photons 
kT, 2) the energy spectral index a (which is related to 
the photon index by the relation T = 1 -I- a), 3) a pa- 
rameter \og{A) related to the Comptonization fraction 
/ (.i.e., the ratio between the number of Compton scat- 
tered photons and the number of seed photons) by the 
relation / = A/{1 + A), and 4) the normalization A^bmc- 

In simple words, the BMC model (which is imple- 
mented in the spectral fitting package Xspec) convolves 
the thermal seed photons and a generic Comptonization 
Green's function producing a power law. As a conse- 
quence, this spectral model generally provides a good 
fit for X-ray spectra of accreting black holes, since their 
continuum is adequately described by a power law. 

The BMC model presents two important advantages 
with respect to the power-law model (PL), which is often 
used to parametrize the Comptonization component: 1) 
Unlike the PL, which is a phenomenological model, the 
BMC parameters are computed in a self-consistent way; 
2) Unlike the PL, the power law produced by BMC does 
not extend to arbitrarily low energies and thus does not 
affect the normalization of the thermal component nor 
the amount of local absorption, which is often present 
around accreting objects. 

2.2. A new scaling method 

Generally, during their transitions from the low-hard 
(LH) to the high-soft (HS) spectral state different GBHs 
show similar spectral variability patterns. Specifically, 
when the X-ray spectra are fitted with the BMC model, 
and r is plotted versus AbmCi the spectral evolution 
appears to be characterized by two saturation levels (the 
lower one corresponding to the LS and the higher one to 
the HS) connected by a monotonically increasing curve. 
Following ST09, this spectral trend can be parameterized 
by the following functional form: 

r(ArBMc) = A-Bln{exp [l - {N^Mc/Nt,)'^] + l} (1) 

Note that ST09 provide in their Equation (10) a slightly 
different functional form that comprises an additional pa- 
rameter D in order to be able to fit both the QPO - F 
and A^BMC " P trends. However, since D = 1 for the pat- 
tern of our interest. Equation (1) coincides with Eq.(lO) 
of ST09. 

By virtue of the similarity of these A'bmc - P trends 
among GBHs, the value of A/bh (and the distance, when 
the A^BMC - P plot is used in combination with the QPO 
- F plot) of any GBH can be obtained from a direct scal- 
ing process. Simply speaking, if Mbh (and the distance) 
is known for a suitable GBH considered as a reference 
system, the black hole mass for any other GBH can be 
determined by horizontally shifting its self-similar func- 
tion until it matches the reference object's pattern (a 
visual explanation of this scaling process is provided in 
Section 3 when this technique is applied to AGNs). 

The physical basis of these scaling techniques is thor- 
oughly explained by ST09 and essentially relies on few 
reasonable considerations that can be summarized as fol- 
lows: (1) the QPO frequency is inversely proportional to 
ATbh (which can be readily understood considering that 
the larger Mbh, the larger the distance from the BH 



of the accretion inflow responsible for the radiation, and 
hence the longer the dynamical timescale of the process) ; 
(2) the BMC normalization is a function of luminosity 
and distance: A'bmc oc L/(P (by definition of the BMC 
model); (3) the luminosity of an accreting black hole sys- 
tem can be expressed as L oc -qM-Qyirn, where r? is the ra- 
diative efficiency and m the accretion rate in Eddington 
units. The self-similarity of the A^bmc - P correlation 
implies that, in the same spectral state, different sources 
have similar values of rj and to. Therefore, the photon 
index is a reliable indicator for the source's spectral state 
independently of its BH mass. 

In simple terms, the necessary steps to derive A/bh 
with this method can be summarized as follows: 

(1) Construct a P — A^bmc plot for a GBH of known mass 
and distance, which will be used as reference (hereafter 
denoted by the subscript r) and for a target of interest 
(denoted by the subscript t). 

(2) Infer the normalization ratio between the target and 
the reference object A^BMC,t/ABMC,r by shifting in the 
F — A^BMC plot the target's pattern until it matches the 
reference one. 

(3) Derive the black hole mass using the following equa- 
tion 

MBH.t = AfBH.r X (A^BMC.t/A^BMC.r) X {dt/drf X fa (2) 

where A/bh.i- is the black hole mass of the reference ob- 
ject, A'^BMCt and A^BMC.r are the respective BMC nor- 
malizations for target and reference objects, dt and dr 
are the corresponding distances, and fa = cos 9r/ cos 9t 
is a geometrical factor that depends on the respective in- 
clination angles and should be included only in the sce- 
nario where the X-ray soft photon emitting region has a 
disk-like geometry. 

Following ST09, we assume GRO J1655-40 as our pri- 
mary reference source, since its system's parameters are 
the most tightly constrained. However, for completeness, 
we will consider all the reference sources and patterns 
used by ST09. The basic properties of these GBHs are 
reported in Table 1. 

In Figure [T] we show the F — A^bmc plots for our three 
reference sources during rise and decay phases of different 
outbursts (the details of these transitio ns are provided in 
ST09 and iTitarchuk fc Shaposhnikovl [20ll . Although 
similar at first sight, the spectral patterns differ from 
source to source as evidenced by the bottom right panel 
of Fig. [TJ However, one must bear in mind that the F — 
A^BMC plot yields the A^BMC,t/ABMC,r ratio, but the A/bh 
estimate also depends on the values of mass and distance 
of the reference source (see Eq. 2). As a consequence, 
reference sources with different patterns (and different 
mass and distance) may lead to similar estimates of A/bh 
(see Section 4.2). 

From Figure [T] we also note that the same source dis- 
plays different spectral patterns during the rise and the 
the decay phase of an outburst. The latter behavior 
simply reflects the observational fact that during any 
outburst every GBH shows an hysteresis loop when the 
spectral properties (in this case the photon index P) are 
plotted versus the its intensity (here parameterized by 
the normalization Abmc)- 

As outlined by ST09, when this method is utilized to 
estimate A^Tbh in a given GBH, one must use the ap- 
propriate reference plot. For example, a F — A^bmc plot 
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Fig. 1. — Photon index measurements plotted versus the BMC normalization values for different transitions of our three reference sources 
GRO J1655-40, GX 339-4, and XTE J1550-564. The thick solid line represents the best fit, whereas the dotted lines are the la confidence 
levels. The small bottom panels describe the data-to-model ratios. Details on the fitting procedure and the values of the functional 
parameters are provided in the text and in Table 2. The right bottom panel shows a comparison between the spectral variability trends of 
the three reference sources on the same scale. 
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obtained during an outburst rising phase for a given tar- 
get source can only be compared to a rising pattern of 
a suitable reference source that shows a similar trend. 
When applying this method to AGNs, given the much 
longer timescales involved, we cannot construct a com- 
plete pattern in the F— iVBMC plot nor determine whether 
the AGN is in a rising or decaying outburst phase. There- 
fore, for completeness, we need to use all the available ref- 
erence patterns and then determine which one provides 
the best match with the values obtained from reverber- 
ation mapping. 

We fitted the X-ray spectral trends with the functional 
form described in Eq uation (1) using the Levenberg- 
Marquardt algorithm (jPress et al.lll997l ). The results of 
the fit, which are reported in Table 2 along with the la 
uncertainties on each parameter, are broadly consistent 
with those of ST09, the main difference being that in our 
fits all parameters are left free to vary. In this way, we 
can determine the uncertainty on each parameter, which 
is then taken into account to estimate the uncertainty on 
the Mbh values. 

A priori, any of the patterns shown in Figure [T] can 
equally well be used to determine Mbh in our sample 
of AGNs. Indeed, all three reference sources have spe- 
cific advantages: GRO J1655-40 has the best determined 
system parameters and therefore reduces the uncertainty 
associated with the estimate of Mbh; GX 339-4 is known 
to be a prototypical GBH from the spectral variability 
point of view with consistent patterns in all the outbursts 
observed; finally, XTE J1550-564 displays the largest 
range of photon indices during its spectral transition, 
allowing to extend the estimate of Mbh to AGNs with 
steep F. However, a visual inspection of Figure [1] reveals 
that the patterns describing the outburst rise phase for 
the three reference sources (top and middle right pan- 
els and bottom left panel) lack the coverage of at least 
one part of the spectral pattern: part of the rising trend 
and all the upper saturation level for GRO J1655, the 
rising trend for GX 339-4, and both saturation levels for 
XTEJ1550-564. On the other hand, the decay phases of 
GRO J1655-40 and GX 339-4 (hereafter GRO J1655 D05 
and GX 339 D03, respectively) are well-sampled across 
the entire range of the spectral pattern, and hence may 
be considered a priori more reliable. Nevertheless, only 
a posteriori, after a comparison of the Mbh estimates 
with the corresponding values obtained from reverbera- 
tion mapping will we be able to assess which reference 
pattern is preferable and in which context. 

3. TEST OF THE SCALING METHOD ON ACTIVE 
GALACTIC NUCLEI 

Despite the large difference in scales, both GBHs and 
AGNs are believed to harbor the same central engine: 
a black hole and an accretion disk/corona that some- 
times produces relativistic jets. There is now mount- 
ing evidence that AGNs may be considered as large- 
scale analogs of GBHs (see, e.g., iKording et al.l 120061 : 
iMcHardv et all l2006t ISobolewska et al.l l2009t ) . There- 
fore, the progress made in the field of GBHs can be in 
principle extended to AGNs (and vice versa). 

Because of their higher brightness (due to their vicin- 
ity) and their shorter variability timescales (direct conse- 
quence of their smaller Mbh), the temporal and spectral 
properties of GBHs are much better known and can be 



used to infer information on their more powerful, extra- 
galactic analogs. In the framework of the AGN-GBH uni- 
fication, it thus appears reasonable to extend to AGNs 
the scaling method described before to determine A/bh- 
To illustrate how this method can be extended to 
AGNs, in Figure [2] we show the F — iVBMC diagram for a 
GBH reference source (thick solid line), and for an hypo- 
thetical AGN with F 1.8 (filled circle). The only basic 
assumption made is that, on a much longer timescale 
than GBHs, the AGN follows a similar spectral pattern 
(light dashed line). The scaling factor necessary to ob- 
tain Mbh, AGN starting from a reference value Mbh.gbh 
is given by the product (iVBMC,t/iVBMC,r) x {dt/d^Y, 
where the first factor is determined by shifting right- ward 
the dashed gray (red, if printed in color) trend until it 
matches the solid darker one. In other words the dif- 
ference in Mbh between the target of interest and the 
reference system is directly related to the amount of the 
shift along the x-axis in the F — A'bmc diagram and to 
the square of the distance ratio. 
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Fig. 2. — A'^BMC - r diagram for a GBH reference source and for 
an hypothetical AGN. The thick continuous Une represents the best 
fitting function of the GBH spectral trend. The arrow illustrates 
how we determine the value of A''bmc t/^BMC r from the value of 
r measured for the AGN. See text for further details. 

3.1. The AGN sample 

We chose a sample of AGNs that fulfills two basic crite- 
ria: (I) the AGNs have BH mass estimates derived from 
a primary method, and (2) possess high quality X-ray 
data, in order to tightly constrain the parameters of the 
BMC model. Nearly 30 object s with Mbh determined 
via reverberation mapping by Pe terson et al.l ([2004) have 
been observed by XMM-Newton, which, thanks to its 
high throughput in the 0.3-10 keV energy band, repre- 
sents the best X-ray satellite for this kind of analysis. 

The main properties of this sample, reported in Ta- 
ble 3, can be summarized as follows: the sample spans 
a redshift range of 0.002-0.234; the values of A^bh en- 
compass nearly 3 orders of magnitude with log(MBH) 
ranging between 6.23 and 9.11 solar masses; the bolo- 
metric luminosities, computed by integrating the spec- 
tra l energy distribution o ver th e 0.001-100 keV interval 
by ^Vasud evan fc FabianI ()2009D , span nearly 5 orders of 
magnitudes ranging between 4 x 10^^ — 1.3 x 10^*" erg s~^ ; 
finally, the Eddington ratio AEdd, obtained from the ra- 
tio between bolometric and Eddington luminosity, covers 
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the range between 0.001-1.14. In summary, not only this 
AGN sample fulfills the 2 basic criteria, but also all the 
relevant physical parameters (with particular emphasis 
for Mbh and the Eddington ratio, which are crucial for 
the comparison between AGN and GBHs) span a con- 
siderably large region of the parameter space, providing 
the ideal framework to test the scaling method. 

3.2. Data reduction 

Since basically all sources are bright and with relatively 
large exposures (the XMM-Newton observation identi- 
fiers and the net exposures are reported in Table 4), we 
restrict our analysis to the EPIC pn data, which provide 
the highest S/N in the energy band of interest. The EPIC 
data were processed in a homogeneous way using the lat- 
est CCD gain values. For the temporal and spectral anal- 
ysis, events corresponding to pattern 0-4 (singles and 
doubles only) were accepted. Arf and rmf files were cre- 
ated with the XMM-Newton Science Analysis Software 
(SAS) 8.0. The recorded events were screened to remove 
known hot pixels and data affected by periods of fiar- 
ing background. In general, the extraction radius used 
for the source spectra and light curves is 30", whereas 
background spectra and light curves were extracted from 
source-free circular regions of 60" extraction radius on 
the same chip as the source. 

The X-ray spectral analysis was performed usi ng the 
XSPEC V. 12.4.0 software package (| ArnaudI IT996I ) . The 
EPIC data have been re-binned in order to contain at 
least 20 counts per channel for the statistic to be valid. 
The errors on spectral parameters are at 90% confidence 
level for one interesting parameter (Ax^ = 2.71). 




model that comprises a BMC model and two absorption 
components, one fixed at the Galactic value and one free 
to vary to mimic the intrinsic local absorption. When 
necessary, one or two Gaussian components were added 
to account for line-like features. 
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Fig. 4. — Distributions of the A/bh values obtained with rever- 
beration mapping (filled histogram), and with the X-ray scaling 
method using as reference GRO J1655-40 D05 (dashed histogram 
with positive slopes) and GX 339-4 D03 (dashed histogram with 
negative slopes). 



We restricted our analysis to the 2-10 keV energy range 
to avoid the complexity related to the soft X-ray band, 
which includes the still debated nature of the soft ex- 
cess and the possible presence of warm absorbers. For 
completeness, we tested whether and how the estimate 
of A/bh is affected by the use of the hard X-ray data only. 
To this end, we first used PKS 0558-504, a bright AGN 
for which we have high quality proprietary XMM-Newton 
data and simultaneous covera ge of the optical/U V band 
with the SWIFT \]WOT fsee lGliozzi etal|[2010l for de- 
tails). We performed a spectral fitting analysis of both 
the extended optical/UV to X-ray range as well as of the 
restricted 2-10 keV energy band. The main difference 
is the resulting temperature (feT ~ 8 — 20 eV when us- 
ing the full range, as opposed to kT ~ 100 eV obtained 
from the 2-10 keV band), which however has a negligi- 
ble impact on the estimate of A/bh • The resulting values 
of A/bh differ by less than 50%. The little impact of kT 
(and of the parameter log (A) related to the Comptoniza- 
tion fraction) on the A^bh estimate was further assessed 
considering the two AGNs of our sample with the high- 
est (Ark 120) and lowest (NCG 3227) Eddington ratio. 
When we arbitrarily varied by 80% the best-fit values of 
kT and log(A), the consequent changes of A/bh were of 
the order of 20-30% or less. We can therefore conclude 
that using the 2-10 keV range does not significantly af- 
fect the determination of A^bh for this sample of AGNs. 
The results of the spectral analysis are summarized in 
Table 4. 



Fig. 3. — Distribution of the Eddington ratio values AEdd = 
^bol/^Edd of the AGNs for which the A/bh was determined with 
the X-ray scaling method. 

4. RESULTS 

4.1. Spectral Analysis 

We carried out a systematic spectral analysis of the 
X-ray data of all AGNs listed in Table 3 using a baseline 



4.2. Estimate of Mbh 

The estimate of A^bh with this scaling method cru- 
cially depends on the properties (A^bh, distance, and 
perhaps inclination angle) of the GBH primary reference 
sources. Although there is a general agreement on the 
GRO J1655-40 system parameters quoted in Table 1, it 
must be noted that the dist ance value has been ques- 
tioned bv lFoellmi et al.l ()2006l ). who proposed a consider- 
ably smaller value {d < 1.7 kpc). However, it was pointed 
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Fig. 5. — Mbh values obtained with the scaling method using different reference sources plotted versus the reverberation mapping values. 
The dashed line is the linear best-fit, the thick solid line indicates the one-to-one correlation; the dotted lines represent the 0.3 dex levels, 
commonly assumed as uncertainty on the reverberation mapping estimates. The geometrical factor Fq is fixed to 1 illustrating a spherical 
geometry scenario. In the bottom-left panels the data points have been shifted along the y axis by a factor of 3. 




Fig. 6. — Mbh values obtained with the scaling method using different reference sources plotted versus the reverberation mapping values. 
The dashed line represents the linear best-fit, the solid line indicates the one-to-one correlation and the dashed lines are the 0.3 dex levels. 
The geometrical factor Fq is determined by the inclination angles reported in Tables 1 and 3 and illustrates a disk-like geometry. 
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out th at the systemic velocity reported bv lFoellmi et alj 
(|2006D is at odds with high er accuracy measurements 
repor ted in the literature ([Gonzalez Hernandez et alJ 
I2OO8D and that the lower value of the distance would 
not allow the companion st ar to fill its Roche lobe 
(jCaballero Garcia et al.l[2007[ ) . In addition, the fact that 
using d = 3.2 kpc ST09 obtain an excellent agreement 
between the values of A/bh determined with the scal- 
ing technique and the corresponding dynamical values 
for several GBHs c asts further doubts on the claim of 
iFoellmi et al.l (|2006f ). In the remainder of the paper, we 
will therefore utilize for the reference sources the system 
parameters provided in Table 1 (if we used d = 1.7 kpc 
for GRO J1655-40 the AGN Mbh values would increase 
by a factor of ■~3.5). 

For 24 sources, the BH mass was estimated using Equa- 
tion (2) and the five patterns of the reference sources 
shown in Fig.Hl Five objects (Mrk 79, NGC 3516, NGC 
4151, PG 1226+023, and PG 1411-^442), perhaps be- 
cause of high intrinsic absorption and/or in refiection- 
dominated state, yielded photon indices too low (F < 
1.4) to allow a direct comparison with the spectral trends 
of the reference sources introduced in Section 2.2. In 
those cases, it was not possible to determine A/bh with 
this new scaling method. The distribution of the loga- 
rithm of the Eddington ratio AEdd = -^boi/^Edd for re- 
maining sample is shown in Figure [3l which illustrates 
that AEdd spans a range of 0.01-1 with the majority of 
the objects populating the 0.03-0.3 interval. 

Since the geometry of the X-ray emitting region is still 
a matter of debate, we have considered two different sce- 
narios: (1) spherical geometry, which is defined by setting 
i^G = 1 in the Mbh equation, and (2) disk-like geometry 
where fa = cos Or/ cos6t and the inclination angles given 
in Tables 1 and 3 are used for the reference sources (for 
GX 339-4 we used 70°) and the AGNs, respectively. 

The distributions of the AGN Mbh estimates obtained 
from the two most well-sampled reference transitions 
(GRO J1655 D05 and GX 339 D03) are shown in Fig- 
ure|4]along with the distribution of values from the rever- 
beration mapping sample. The apparent general agree- 
ment between the distributions is formally confirmed 
by a Kolmogorov-Smirnov (KS) test, which indicates 
that the distributions obtained with this new X-ray scal- 
ing method are indistinguishable from the reverberation 
mapping one (for both reference sources the KS statistic 
is ~0.19 with an associated probability of ^81%). We 
also report the individual values of Mbh in Table 5 to 
allow a direct quantitative comparison with the corre- 
sponding values from reverberation mapping. The er- 
rors quoted are computed adding in quadrature the un- 
certainties on the A^BH and distance of the reference 
sources (given in Table 1) and the uncertainty on the 
(-^BMC,t/A^BMC,r), which in turn depends on the 90% 
confidence errors of the X-ray spectral analysis on both 
A'bmc and F (in Table 4) as well as on the uncertain- 
ties of the functional parameters describing the reference 
patterns in the F — iVBMC diagram. 

The results in Table 5 suggest that this method pro- 
vides A^BH estimates that are well constrained (with an 
average percentage error of '~45% when using the GRO 
J1655 D05 transition and ~85% for GX 339 D03) and 
in general agreement with the corresponding values de- 
termined via reverberation mapping (the average ratio 



between the Mbh values determined with this method 
and those from reverberation mapping is 1.5 ± 0.4 for 
GRO J1655 D05 and 2.2 ± 0.5 for GX 339 D03). 

Reasonable values for AGN Mbh are derived using all 
the reference spectral transitions available, as demon- 
strated by the values reported in Table 6. The ith col- 
umn of this table can be considered as the average depar- 
ture from the mass ratio of 1 and hence as a proxy of the 
discrepancy between the values of AfBH determined with 
X-ray scaling method and the corresponding reverbera- 
tion mapping ones. More specifically, (max(AfBH, ratio)) 
is obtained by inverting the Mbh ratios that are lower 
than the unity (as an example, for Fairall 9, instead of 
using 0.3 as reported in columns 3 of Table 5 we use 
3.3=1/0.3) and then averaging over all Mbh ratios. In 
this way, values much lower than unity will not cancel 
out with values much larger than unity. The 3'"'' col- 
umn of Table 6 reports the regular average ratio be- 
tween AfBH values determined with the two methods. 
This provides information about the tendency of differ- 
ent reference patterns to either overestimate or underes- 
timate A^BH with respect to the reverberation mapping 
values. For example, the spectral patterns corresponding 
to a decay phase of the outburst tend to slightly overesti- 
mate A/bh when a spherical geometry is assumed and to 
slightly underestimate it when a disk geometry is chosen. 
On the other hand, the rising patterns have a tendency 
to underestimate A^bh by a factor of ~2-3 when a spher- 
ical geometry is considered and by a larger amount when 
a disk-like geometry is used. 

We can conclude that all the reference patterns used 
with a quasi-spherical illuminating soft photon emission 
provide reliable estimates of Mbh in our AGN sample, 
with GRO J1655 DOS and XTE J1550 R98 showing the 
lowest and highest average discrepancy with respect to 
the reverberation mapping estimates {'^2 and ~6, re- 
spectively). 

4.3. Correlation Analysis 

An additional direct way to compare the A^bh es- 
timates from this method with the corresponding val- 
ues determined via reverberation mapping is to plot one 
quantity versus the other and perform linear correlation 
tests. The A^bh estimates derived from all the refer- 
ence transitions are plotted against the corresponding 
values obtained from the reverberation mapping in Fig- 
ure^ (spherical geometry) and Figure |6] (disk- like geome- 
try), where the dashed line represents the linear best-fit, 
the solid line indicates the one-to-one correlation and the 
dotted lines are the 0.3 dex levels, commonly assumed 
as uncertainty on the reverberation mapping estimates. 
The vertical error-bars account for the uncertainties in 
the BH mass and distance of the reference sources, for 
the uncertainty of the functional reference and spectral 
parameters, as explained in the previous section. 

A visual inspection of these figures reveals a general 
agreement between the Mbh values obtained with these 
two completely independent techniques. This consis- 
tency is more evident in the spherical geometry scenario 
shown in Figure [5l In the bottom right panel of this 
figure, we show that an excellent agreement can be ob- 
tained also using the XTE J1550-564 transition as a 
reference, provided that the Mbh values are systemat- 
ically increased by a factor of 3. The tight correlation 
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observed in Figures [5] and [6] is formally confirmed by 
non-parametric correlation analyses based on the Spear- 
man's p and Kendalls's r rank correlation coefficients as 
well as by a linear correlation fit using the routine f itexy 
([Press et al.lll997t ) that accounts for the errors not only 
on the y-axis but along the x-axis as well. 

In Table 7 we report the results of this analysis. From 
the results of the Spearman's and Kendalls's analyses we 
infer that the Mbh values obtained with these two dif- 
ferent method are always correlated at high significance 
level regardless of the reference transition used. This 
conclusion is further supported by the fact that in all 
cases the best-fit line's slope is either consistent or very 
close to unity (and inconsistent with the zero value at 
more than lOcr level) and the intercept is close to zero. 

5. DISCUSSION AND CONCLUSIONS 

In this work, we have explored whether a new X- 
ray scaling method recently introduced to determine the 
mass of stellar black holes can be successfully extended to 
constrain the mass of supermassive black holes in AGNs. 
To this aim we have utilized a sample of 24 AGNs with 
r = 1.56 - 2.11 and AEdd = 0.01 - 1, and with Mbh pre- 
viously determined via reverberation mapping and with 
good-quality archival XMM-Newton data. The main re- 
sults can be summarized as follows. 

• This novel method, which is based on the spectral 
fitting of X-ray data with the BMC Comptoniza- 
tion model and on the similarity of the X-ray spec- 
tral behavior in BHs regardless of of their mass, 
appears to be a robust estimator of Mbh in AGNs. 

• The values of Mbh determined with this method 
are well constrained and in good agreement (within 
a factor of ~2-3) with the corresponding values ob- 
tained with the reverberation mapping technique. 

• For AGNs with F comprised in the 1.4-1.95 range 
the most reliable reference appears to be the out- 
burst decay pattern of GRO J1655-40 in combi- 
nation with spherical scenario (i.e., with Fq — 1 
in Eq. 2). For F < 2.1 both decay and rise pat- 
terns of GX 339-4 and spherical geometry provide 
reliable estimates of Mbh- Finally, for steeper pho- 
ton indices up to F '-^ 3 only the rising pattern of 
XTE J 1550-564 can be used, keeping in mind its 
tendency to underestimate Mbh by a factor of ^3. 

• The A/bh values appear to be systematically under- 
estimated by a factor of 2-3 when a disk-like geome- 
try is used instead of a quasi-spherical geometry for 
the soft photon supply. This result lends support to 
the hypothesis that the soft photon emission area 
is quasi-spherical. 

In summary, our analysis reveals that this X-ray scal- 
ing technique yields black hole masses that are in full 
agreement with the accepted values obtained from the 
reverberation mapping. This is consistent with the hy- 
pothesis that the X-ray spectrum produced by the Comp- 
tonization process is self-similar and its shape is indepen- 
dent of Mbh, and lends further support to grand unifi- 
cation model between AGNs and GBHs. 



Before concluding it is important to briefly discuss the 
range of applicability of this technique and its main dif- 
ferences with respect to most commonly used methods 
for black hole mass determination. First, this is one 
of the few truly scale-independent techniques to esti- 
mate Mbh- The other two scale-independent methods 
are the one based on the so called "fundamental plane" 
of b lack holes (BHs) int roduced bv lMerloni et all ()2003f ) 
and iFalcke et al.l ([2004'), where Mbh is related to both 
the X-ray and radio luminosities, and the one based on 
the scaling of te mporal breaks in X-ray power spectral 
density analyses (jMcHardv et al.ll2006i ). However, the 
fundamental plane is affected by large scattering and 
is limited to BH systems characterized by radio emis- 
sion, whereas the method based on X-ray timing requires 
high-quality evenly-sampled X-ray light curves spanning 
a baseline of more than a decade and therefore is neces- 
sarily limited to very few selected AGNs. On the other 
hand, this X-ray scaling novel technique only needs mod- 
erately good X-ray spectra, which are now available for a 
very large number of AGNs and BH systems in general. 
Second, unlike most indirect methods, this X-ray scaling 
method is completely independent of any assumption on 
the BLR nature/geometry or host galaxy characteristics, 
and can therefore complement several optically-based 
studies by providing independent estimates of A/bh, and 
expand them when the optical properties are uncertain 
or unavailable. 

In principle, this technique can be applied to any ac- 
creting black hole with bulk X-ray emission produced via 
Comptonization and for which the photon index falls in 
a range covered by the GBH reference patterns in the 
F — A^BMC diagram. However, we caution that a blind 
application of this technique may lead to erroneous re- 
sults. This will occur whenever the genuine properties of 
the primary Comptonized X-ray emission are not prop- 
erly determined, for example because effects associated 
with strong absorption, reflection, or emission from an 
additional component such as a relativistic jet are not 
carefully accounted for. Particular care should be used 
when dealing with low accreting systems, because this 
method has not been tested below Asdd = 0.01 and be- 
cause at very low values of m a radiatively inefficient 
accretion mode (possibly accompanied by outflows and 
jets) is likely to take place leading to the bulk of the 
X-rays to be domin ated by jet emission as suggested by 
lYuan fc Cuil ([200l . 

It is also worth mentioning that AGNs are generally 
characterized by high flux variability, which is often ac- 
companied by spectral variability. The vast majority 
of bright radio-quiet AGNs show steeper spectra when 
they are brighter and hence their evolution is in quali- 
tative agreement with the rising trend in the F — A^bmc 
plot, providing strong support for the physical basis of 
our methodology. Additionally, AGN flux and spectral 
changes are generally not so extreme and therefore the 
value obtained with different observations will be within 
the uncertainty of the method. Finally, occurrences 
of extreme spectral changes are usually associated with 
changing in the local obscuration component (for exam- 
ple from Compton-thick to Compton-thin scenario) and 
therefore do not reflect the intrinsic change of the Comp- 
tonized primary component. A possible notable excep- 
tion is NGC 4051, which, in virtue of its small SMBH, 
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TABLE 1 

Characteristics of reference sources 



Name 


Mbh 


d 






(Mq) 


(kpc) 


(deg) 


(1) 


(2) 


(3) 


(4) 


GRO J1655-40 


6.3 ±0.3 


3.2 ±0.2 


70 ± 1 


GX 339-4 


12.3 ± 1.4 


5.7 ±0.8 


> 45 


XTE J1550-564 


10.7 ± 1.5 


3.3 ±0.5 


72 ±5 



Columns Table 1: 1= Source na me. 2= Black hol e ma ss. 3= Distance. 4= Inclin ation angle. Notes: For GRO J1655-40 the quoted 
values of A/bh, d, and i are from G reene et al.l l|200ll ) and lHiellming fc Rupenl HMEb. Fo r both GX 339-4 and XTE J1550-56 4 Mbh and 
d are from ST09, whereas the values of i are respectively from from lKolehmainen Sz Pond l|201(3 'l and from lOrosz et aP l|2002l ). 



shows extreme flux and spectral changes that might re- 
flect genuine changes of the primary X-ray component. 
This is indeed the only source for which a complete trend 
in the T — Nbmc plot can be constructed. Preliminary 
results based on numerous XMM-Newton and Chandra 
observations suggest that NGC 4051 follow a trend con- 
sistent with GBH references and that the scaling method 
successfully determine its mass (Chekhtman & Titarchuk 
2011, in preparation). 

In the near future, we plan to systematically apply 
this X-ray scaling method to different categories of ob- 
jects, whose Mbh values are currently debated, such as 



ultraluminous X-ray sources, narrow line Seyfert 1 galax- 
ies, non-jet-dominated radio loud AGNs, and powerful 
quasars to test the upper end of the Mbh ~ cr^, relation- 
ship. 
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TABLE 2 

Parametrization of r — ATbmc reference patterns 



Transition 


A 




B 






Xl (dof) 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


GRO J1655 DOS 


1.96 ±0.02 


0.42 


±0.02 


0.023 ± 0.001 


1.8 ±0.2 


1.5 (20) 


GRO J1655 R05 


2.35 ± 0.04 


0.74 


±0.04 


0.131 ±0.001 


1.0 ±0.1 


1.8 (34) 


GX339 DOS 


2.13 ± 0.03 


0.50 


± 0.04 


0.0130 ± 0.0002 


1.5 ±0.3 


0.9 (40) 


GX339 R04 


2.10 ±0.03 


0.46 


±0.01 


0.037 ±0.001 


8.0 ± 1.5 


2.3 (24) 


XTE J1550 R98 


2.96 ± 0.02 


2.8 


±0.2 


0.055 ± 0.010 


0.4 ±0.1 


2.7 (49) 



Columns Table 2: 1= Reference source spectral transition. 2= Parameter of the functional form described in Eq.(l) that is responsible 
for the rigid translation of the spectral pattern along the y-axis. 3= Parameter characterizing the lower saturation level of the pattern. 
4= Parameter responsible for the translation of the spectral pattern along the x-axis. 5= Slope of the spectral pattern. 6= Reduced 
and degrees of freedom. Notes: The quoted parameter errors are the lir uncertainties obtained from fitting the spectral data using the 
Levenberg-Marquardt algorithm. 
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TABLE 3 
Sample properties 



Name z (i) log(MBH) log(ibol) ^Edd 

(dcg) (Mo) (erg s-i) (Lbol/^Edd) 



(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


3C 120 


0, 


,0330 




7.74 


45.3 


0.305 







.0561 


48+24 


8.46 


45.2 


0.047 


Ark 12n 





0327 


99+16 


8.18 


45.3 


0.111 




n 
u 






8 41 


44 8 




Mrk 110 





.0353 


1 1 +9 


7.40 


45.1 


0.433 


TV/TfL- 97Q 


n 
U 


.uouo 






^o.u 


910 
U.ZiU 


Mrk 335 





.0258 


o_4 


7.15 


45.3 


1.130 


Mrk 509 





.0344 


5+* 

'-'-2 


8.16 


45.2 


0.095 


Mrk ^QO 

1V±1 JV \Jij\J 





.0264 


10+9 


7.68 


43.8 


0.010 


Mrk 79 





.0222 


01+23 
•^^-13 


7.72 


44.3 


0.031 


NOC 'H227 





.0039 


71+17 
' -'--24 


6.88 


42.9 


0.001 


NGC 3516 


0, 


.0088 


30 


7.50 


43.5 


0.006 


NGC 3783 


0, 


.0097 




7.47 


44.2 


0.043 


NGC 4051 


0, 


.0023 




6.23 


42.6 


0.016 


NGC 4151 


0, 


.0033 




7.12 


44.0 


0.056 


NGC 4593 


0, 


.0090 


30 


6.99 


43.7 


0.037 


NGC 5548 


0, 


.0172 


45lf8 


7.65 


44.3 


0.024 


NGC 7469 


0, 


.0163 




7.09 


44.8 


0.369 


PG 0052+251 


0, 


.1550 


5t* 


8.57 


45.8 


0.148 


PG 0844+349 


0, 


.0640 


8i| 


7.97 


45.4 


0.233 


PG 0953+414 


0, 


.2341 


9+2 
^-1 


8.44 


46.5 


0.892 


PG 1211+143 


0, 


.0809 




8.16 


45.7 


0.260 


PG 1226+023 


0, 


.1583 




8.95 


47.1 


1.140 


PG 1229+204 


0, 


.0630 


12tf 


7.86 


44.9 


0.082 


PG 1307+085 


0, 


.1550 


5+5 
3-2 


8.64 


45.6 


0.066 


PG 1411+442 


0, 


.0896 


3t? 


8.65 


45.4 


0.041 


PG 1426+015 


0, 


.0865 


13t^^ 


9.11 


45.6 


0.024 


PG 1613+658 


0, 


.1290 


15^12 


8.45 


45.9 


0.221 


PG 2130+099 


0, 


.0630 


6+* 
"-3 


7.58 


45.0 


0.018 



Columns Table 3: 1= Source name. 2= Redshift. 3= Inclination angle. 4= Black hole mass determined via reverberation mapping. 
5= Bolomctric luminosity. 6= Eddington rati o. Notes: The AGN distances were computed from z assuming Hq = 71 km s~^ Mpc"!, 
f^A = 0.73, and Vim = 0.27 llSpergel et al.|200 3l). For nearby AGNs we used the redshift-independent distances provided by the NASA/IPAC 
extragalactic database (NED). The quoted values of the inclination angle are from Bian & Zhao (2002); for objects without estimate 
of inclination angle (Mrk 279, NGC 3516, NGC 4593) we assumed a value of 30° typical for Seyfert galaxies. M bh values are from 

ticrson ct al. (2004) with the exception of NGC 3227, NG C 3516, NGC 4051. NGC 5548 l|De'nnev et al.ll20Tol V and PG 2130+099 
'ier et al.ii2008i) . Values for (ibol) and ^Edd from 1 Vasudevan &c FabianI (120091) . 
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TABLE 4 
Spectra Results 



Name 

(1) 


XMM obsid 

(2) 


Net exposure 

(ks) 
(3) 


kT 

(keV) 
(4) 


Log{A) 
(5) 


r 

(6) 


-^BMC 
(10-4) 

(7) 


Xr(dof) 

(8) 


3C 120 


0152840101 


79.1 


0.4 ±0.1 


0.7 ±0.1 


l.D9_Q 02 


4 5+01 


1.02(1211) 


3C 390.3 


0203720201 


35.0 


0.1 ±0.1 


0.7 




9 i;+0.2 
^•^-0.2 


1.13(817) 


Ark 120 


0147190101 


55.9 


0.4 ±0.4 


0.5 ±0.1 




4 2+01 


1.07(1000) 


Fairall 9 


0101040201 


26.0 


0.1 ±0.1 


0.5 ±0.1 


1.80«;;^ 


1 0+01 


1.05(866) 


Mrk 110 


0201130501 


32.9 


0.3 ±0.1 


0.6 ±0.3 


1.68lHi 


9 7+1.0 
^■'-0.1 


0.99(1290) 


Mrk 279 


0302480501 


34.9 


0.05 


0.95 


1.801°;°! 


1 5+0-1 


1.05(1252) 


Mrk 335 


0306870101 


83.3 


0.3 ±0.1 


0.4 ±0.1 


1 qs+O-03 
^•^^"-0.03 


9 9+0.2 
^■^-0.1 


0.99(1364) 


Mrk 509 


0306090401 


44.9 


0.3 ±0.1 


0.6 ±0.1 


1 fiQ+002 
i.ua_Q Q2 


^•^-0.1 


1.02(1468) 


Mrk 590 


0201020201 


39.5 


0.4 ±0.1 


0.8 ±0.2 


1 i;4+00B 


o-6iri 


1.03(838) 


Mrk 79 


0502091001 


49.0 


0.11 


0.7 


1 99+0.02 
J^-^^_0.02 




1.02(1065) 


NGC 3227 


0400270101 


96.1 


0.08 


0.5 






0.98(1574) 


NGC 3516 


0107460701 


81.8 


0.2 ±0.1 


0.5 


-, 09+0.02 
^•'-'^-0.01 


-1 Q+O.l 

1-^-0.1 


1.06(1411) 


NGC 3783 


0112210501 


89.5 


0.2 ±0.1 


0.8 


1 56+001 


4 2+°-^ 


1.19(1325) 


NGC 4051 


0606320101 


45.7 


0.1 ±0.1 


0.1 ±0.1 


1 79+0.004 
-^•'"'-0.004 


1 5+0-1 


0.99(1059) 


NGC 4151 


0143500301 


12.4 


0.01 


2.0 


1 Qc+O.Ol 

1.00_g Q]^ 


o 4+0.1 
•^■^-0.1 


1.30(1593) 


NGC 4593 


0109970101 


5.0 


0.4 ±0.1 


0.6 ±0.2 


1 cc+O.OS 
^•°>'-0.09 


3 6+°-4 


0.94(648) 


NGC 5548 


0089960301 


50.5 


0.10 


0.6 ±0.1 




9 c+0.1 
^•^-0.1 


1.01(1494) 


NGC 7469 


0207090101 


58.3 


0.2 ±0.1 


1.2 


1.861°;°! 


9 0+0.1 


1.15(902) 


PG 0052+251 


0301450401 


8.2 


0.3 ±0.2 


0.5 


1 7c:+0.08 


0.6l°ii 


0.93(243) 


PG 0844+349 


0103660201 


12.7 


0.3 ±0.1 


0.4 ±0.4 


9 nn+0.12 


6+°° 
"■"-0.1 


0.96(297) 


PG 0953+414 


0111290201 


11.2 


0.1 ±0.1 


0.1 


2 00+°-°'^ 


5+°-l 
"-"-0.1 


0.74(170) 


PG 1211+143 


0502050101 


45.6 


0.1 ±0.1 


2.0 


9 1-1+0.03 
-0.03 


0.41°;^ 


0.90(659) 


PG 1226+023 


0414190101 


45.7 


0.10 


0.5 


, on+O.Ol 
l.>39_() 


6.4tg;l 


1.04(1595) 


PG 1229+204 


0301450201 


16.9 


0.05 


1.0 


1 01+0.06 

-^•^-^-0.06 


3+°-l 
"-■^-0.1 


0.94(253) 


PG 1307+085 


0110950401 


10.6 


0.10 


0.6 


1 60+0-10 


l+o-i 
"-I-0.1 


0.91(98) 


PG 1411+442 


0103660101 


22.4 


0.07 


2.0 


1 00+1-°° 


3.41SJ 


0.92(41) 


PG 1426+015 


0102040501 


5.4 


0.03 


1.2 


1 97+0.06 
-0.06 


o.4iS;l 


0.88(214) 


PG 1613+658 


0102040601 


5.0 


0.04 


1.2 


1 Q4+0-29 
i-y*_0.26 


Q + 0.3 


1.05(14) 


PG 2130+099 


0150470701 


21.9 


0.06 


0.7 


1 70+0.05 

^•'•^-0.06 


2+°i 


0.98(339) 



Columns Table 4: 1= Source name. 2= XMM-Newton observation identifier. 3= EPIC pn net exposure time in ks. 4= Temperature 
of thermal photon source in kcV. 5= Logarithm of the A parameter, which is related to the Comptonization fraction / by the relation 
/ = A/{1 + A). 6= 2-10 keV photon index. 7= Normalization of the Comptonization model in units of (L/IO^® erg s^i)(10 kpc/a!)-^. 8= 
Reduced degrees of freedom. Notes: In the cases, where the error procedure in XSPEC was not able to provide the 90% confidence 

uncertainty for kT or log(yl) due to low S/N data or spectral degeneracy, we fixed the parameters at their best-fit value. The upper limit 
of log(j4) was fixed at the value of 2. 
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TABLE 5 
Black hole mass estimates 



Name GRO J1655 DOS GX 339 D03 

Sphere Disk Sphere Disk 

logMBH scal/RM logMee scal/RM logMee scal/RM logMeH scal/RM 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



3C 120 


8, 


,36 


± 


0, 


,17 


4, 


,1 


7, 


,91 


± 


0, 


,17 


1 


.4 


8.53 


± 


0, 


,21 


6.1 


8, 


,08 


± 


0, 


,21 


2, 


,2 


3C 390.3 


8. 


,61 


± 


0, 


,20 


1, 


.4 


8, 


,32 


± 


0, 


,20 





.7 


8.80 


± 


0, 


,27 


2.2 


8, 


,50 


± 


0, 


,27 


1, 


.1 


Ark 120 


8, 


,17 


± 


0, 


,16 


1, 


,0 


7, 


74 


± 


0, 


16 


0, 


.4 


8.32 


± 


0, 


,17 


1.4 


7, 


,89 


± 


0, 


,17 


0, 


,5 


Fairall 9 


7, 


,91 


± 


0, 


,24 


0, 


,3 


7, 


,46 


± 


0, 


,24 


0, 


.1 


8.06 


± 


0, 


,27 


0.5 


7, 


,61 


± 


0, 


,27 


0, 


,2 


Mrk 110 


8. 


,24 


± 


0. 


,29 


7, 


,0 


7. 


,78 


± 


0. 


,29 


2, 


.4 


8.42 


± 


0, 


,28 


10.5 


7, 


,96 


± 


0, 


,28 


3, 


.7 


Mrk 279 


7. 


,71 


± 


0. 


,25 


1, 


,5 


7. 


,32 


± 


0. 


,25 


0, 


.6 


7.87 


± 


0, 


,17 


2.1 


7, 


,47 


± 


0, 


,17 


0, 


,8 


Mrk 335 






























7.65 


± 


0, 


,18 


3.1 


7, 


,18 


± 


0, 


,18 


1, 


,1 


Mrk 509 


8, 


,34 


± 


0, 


17 


1, 


,5 


7, 


,89 


+ 


0, 


17 


0, 


,5 


8.51 


+ 


0, 


,21 


2.3 


8, 


,05 


± 


0, 


,21 


0, 


,8 


Mrk 590 


7. 


,51 


± 


0. 


,34 


0, 


.7 


7, 


,04 




0, 


,34 


0, 


.2 


7.82 


± 


0, 


,40 


1.4 


7, 


,37 


± 


0, 


,40 


0, 


,5 


NGC 3227 


6. 


,69 


± 


0, 


,17 


0, 


,6 


6, 


,24 


± 


0, 


,17 





.2 


6.97 




0, 


,20 


1.2 


6, 


,99 


± 


0, 


,20 


1, 


,3 


NGC 3783 


7, 


,52 


± 


0, 


19 


1, 


,1 


7, 


,55 


± 


0, 


19 


1, 


,2 


7.79 


± 


0, 


,73 


2.1 


7, 


,50 


± 


0, 


,73 


1, 


,1 


NGC 4051 


5. 


,77 


± 


0, 


,15 


0, 


,3 


5, 


,48 




0, 


,15 


0, 


.2 


5.94 


± 


0, 


,17 


0.5 


5, 


,63 


± 


0, 


,17 


0, 


,2 


NGC 4593 


7. 


,28 


± 


0, 


,28 


2, 


,0 


6, 


,98 


± 


0, 


,28 


1, 


.0 


7.47 


± 


0, 


,31 


3.0 


7, 


,07 


± 


0, 


,31 


1, 


,2 


NGC 5548 


7, 


,62 


± 


0, 


16 


0, 


,9 


7, 


,21 


± 


0, 


16 


0, 


,4 


7.81 


± 


0, 


,22 


1.5 


7, 


,49 


± 


0, 


,22 


0, 


,7 


NGC 7469 


7, 


,23 


± 


0, 


,16 


1, 


,4 


6, 


,91 




0, 


,16 


0, 


.7 


7.39 


± 


0, 


,16 


2.0 


6, 


,94 


± 


0, 


,16 


0, 


,7 


PG 0052+251 


8. 


,85 


± 


0. 


,23 


1, 


,9 


8. 


,41 


± 


0. 


,23 


0, 


.7 


9.01 


± 


0, 


,31 


2.8 


8, 


,54 


± 


0, 


,31 


1, 


,0 


PG 0844+349 






























8.04 


+ 


0, 


,20 


1.2 


7, 


,58 


± 


0, 


,20 


0, 


,4 


PG 0953+414 






























9.07 


± 


0, 


,25 


4.3 


8, 


,61 


± 


0, 


,25 


1, 


,5 


PG 1211+143 






























7.83 


± 


0, 


,22 


0.5 


7, 


,36 


± 


0, 


,22 


0, 


,2 


PG 1229+204 


7. 


,47 


± 


0. 


,20 


0, 


,4 


7. 


,01 


± 


0. 


,20 


0, 


.1 


7.64 


± 


0, 


,22 


0.6 


7, 


,18 


± 


0, 


,22 


0, 


,2 


PG 1307+085 


8. 


,32 


± 


0, 


,35 


0, 


,5 


7, 


,85 


± 


0, 


,35 





.2 


8.54 


± 


0, 


,42 


0.8 


8, 


,08 


± 


0, 


,42 


0, 


,3 


PG 1426+015 






























8.31 


± 


0, 


,19 


0.2 


7, 


,85 


± 


0, 


,19 


0, 


,1 


PG 1613+658 


8, 


13 


± 


0, 


,18 


0, 


,5 


7, 


66 


+ 


0, 


,18 





,2 


8.34 


+ 


0, 


,41 


0.8 


7, 


,89 


± 


0, 


,41 


0, 


,3 


PG 2i:-!()+09!) 




,()() 




0, 


,20 


1 


,0 






± 


0, 


,20 


0, 


, i 


7.7() 


+ 


0, 


,21 






,29 




0, 


,21 


0, 


,5 



Columns Table 5: 1= Source Name. 2= Mbh in solar units computed using the GRO J1655 D05 transition in the case of spherical 
geometry. 3= Ratio between the Mbh in column 2 and the corresponding value determined via reverberation mapping. 4= Same as 
column 2 but for a disk-like geometry. 5= Ratio between the Mbh in column 3 and the corresponding value determined via reverberation 
mapping. 6= Mbh in solar units computed using the GX 339 DOS transition in the case of spherical geometry. 7= Ratio between the Mbh 
in column 6 and the corresponding value determined via reverberation mapping. 8= Same as column 6 but for a disk-like geometry. 9= 
Ratio between the Mbh in column 8 and the corresponding value determined via reverberation mapping. Notes: For five sources (Mrk 
335, PG 0844+349, PG 0953+414, PG 1211+143, PG 1426+015), the photon index is too steep (F > 1.96) to determine Mbh using the 
decay phase of GRO J1655-40. 
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TABLE 6 
Statistical analysis results 



Transition 

(1) 


Geometry 


(MBH,scal/MBH.RM> ( 


maa;(MBH, ratio)) 
(4) 


GRO DOS 


Sphere 


1.5 ±0.4 


2.2 ±0.3 




Disk 


0.6 ±0.1 


3.3 ±0.5 


GRO R05 


Sphere 


0.4 ±0.1 


5.7 ±1.4 




Disk 


0.16 ±0.03 


14.9 ±4.1 


GX D03 


Sphere 


2.2 ±0.5 


2.7 ±0.4 




Disk 


0.8 ±0.2 


3.1 ±0.7 


GX R04 


Sphere 


0.7 ±0.1 


3.3 ±0.6 




Disk 


0.3 ±0.1 


7.7 ± 1.7 


XTE R98 


Sphere 


0.4 ±0.1 


5.8 ± 1.2 




Sphere-shift 


1.1 ±0.2 


2.5 ±0.4 




Disk 


0.13 ±0.02 


16.6 ±3.8 



Columns Table 6: 1= Reference source spectral transition. 2= Geometry of the soft photon emission region. 3= Ratio of the Mbh 
values determined with the scaling method over the corresponding values from the reverberation mapping averaged over the entire sample. 
4= Maximum average ratio of the Mbh values which provides an estimate of the average discrepancy (see Section 4.2 for details). Notes: 
The "sphere-shift" case corresponds to the spherical corona scenario where the Mbh values obtained from the scaling technique have been 
multiplied by a factor of 3. The large values of (maa; (Mbh, ratio)) for patterns corresponding to rising phase outbursts (especially in the 
cases of disk geometry) reflect the systematic underestimate of Mbh not the true departure from the black hole mass ratio of 1. 



TABLE 7 
Correlation analysis results 



Transition Geom. Slope Intercept Spearman (Prob.) Kendall (Prob.) 



(1) 


(2) 


(3) 


(4) 




CS) 






CS) 




GRO D05 


Sphere 


1.09 ± 0.07 


-0.6 ± 0.6 


0.78 


(9 X 10- 


5-1 


0.59 


(4 X 10- 






Disk 


1.03 ± 0.06 


-0.6 ± 0.5 


0.77 


(1 X 10- 


4-1 


0.57 


(5 X 10- 




GRO R05 


Sphere 


0.95 ±0.07 


-0.1 ±0.5 


0.78 


(5 X 10" 


6\ 


0.62 


(2 X 10- 






Disk 


0.89 ± 0.07 


-0.1 ±0.5 


0.77 


(1 X 10- 


5\ 


0.59 


(6 X 10- 




GX D03 


Sphere 


0.93 ± 0.08 


0.7 ±0.6 


0.79 


(4 X 10- 


6\ 


0.63 


(2 X 10- 






Disk 


0.85 ± 0.07 


0.9 ±0.6 


0.76 


(1 X 10- 


5\ 


0.58 


(7 X 10- 




GX R04 


Sphere 


1.03 ±0.07 


-0.5 ± 0.5 


0.84 


(5 X 10- 


7\ 


0.65 


(1 X 10- 






Disk 


0.94 ± 0.06 


-0.2 ±0.5 


0.82 


(2 X 10- 


6\ 


0.64 


(2 X 10- 




XTE R98 


Sphere 


0.97 ±0.09 


-0.4 ±0.7 


0.80 


(2 X 10- 


6\ 


0.62 


(2 X 10- 


') 




Sphere-shift 


0.99 ± 0.08 


-0.03 ±0.6 
















Disk 


0.88 ± 0.09 


-0.1 ±0.7 


0.80 (2 X IO- 




0.62 (2 X IO- 





Columns Table 7; 1= Reference source spectral transition. 2= Geometry of the soft photon emission region. 3= Best-fit slope obtained 
fitting the data in Figures 3 and 4 with a straight line. 4= Best-fit intercept. 5= Spearman's p rank correlation coefficient and related 
chance probability. 6= Kcndalls's r rank correlation coefficient and related chance probability. Notes; The quoted parameter errors for 
slope and intercept account for the uncertainties on both x and y. The uncertainty on the average mass ratio represents the spread around 
the mean a/n, where n is the number of AGNs. For the "sphere-shift" case the values of the Spearman and Kendall analyses have been 
omitted since they coincide with those of the "sphere" case. 



